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We propose that “standard pings”, brief broadband radio impulses, can be used to study the three-dimensional 
clustering of matter in the Universe even in the absence of redshift information. The dispersion of radio waves 
as they travel through the intervening plasma can, like redshift, be used as a cosmological distance measure. 
Because of inhomogeneities in the electron density along the line of sight, dispersion is an imperfect proxy for 
radial distance and we show that this leads to calculable dispersion-space distortions in the apparent clustering of 
sources. Fast radio bursts (FRBs) are a new class of radio transients that are the prototypical standard ping and, 
due to their high observed dispersion, have been interpreted as originating at cosmological distances. The rate of 
fast radio bursts has been estimated to be several thousand over the whole sky per day and, if cosmological, the 
sources of these events should trace the large-scale structure of the Universe. We calculate the dispersion-space 
power spectra for a simple model where electrons and FRBs are biased tracers of the large-scale structure of 
the Universe and we show that the clustering signal could be measured using as few as 10 000 events. Such a 
survey is in line with what may be achieved with upcoming wide-held radio telescopes. 


Introduction .—The clustering of matter on large scales has 
been heralded as the next great probe of the Universe after the 
cosmic microwave background (CMB). The large-scale struc¬ 
ture, in principle, contains far more information than the CMB 
because it can be studied in three dimensions. Traditionally, 
the redshift of spectral lines of “standard atoms”, caused by 
the Hubble flow, has been employed to measure radial dis¬ 
tance and provide the third dimension. Such redshift surveys 
have now aggregated the positions of millions of galaxies into 
three dimensional density maps, exposing a rich structure of 
clusters, filaments, walls, and voids. 

However, redshift is not the only measure of radial distance. 
Standard candles and standard rulers can be used to estimate 
distance using brightness and angular size, respectively. In 
the same way a “standard ping”, a brief broadband radio im¬ 
pulse, can be used to measure dispersion through the char¬ 
acteristic wavelength-squared delay of electromagnetic sig¬ 
nals. The dispersion is proportional to the amount of plasma 
along the line of sight and can thus be used to estimate ra¬ 
dial distance. This technique, routinely used for pulsars in the 
Galaxy, can be used to define a dispersion distance to cosmo¬ 
logical sources. 

Dispersion is readily measured from fast radio bursts 
(FRBs), a new class of millisecond duration radio transients 
first reported by Lorimer et al. Ill], which are the archetype 
of a standard ping. Ten such bursts have been reported to 
date 12-0] . the total event rate has been estimated to be 

3.3^2 5 X 10^ per sky per day Qj. Dispersion is characterized 
by a line-of-site integral of the free electron density known the 
dispersion measure (DM). FRB events are observed to have 
dispersion measures of order ^lOOOpc/cm^, greatly in ex¬ 
cess of the galactic expectation of order ~100pc/cm^ (de¬ 
pending on galactic latitude) Qj. They thus appear to be of 
extragalactic origin, with their dispersion produced either by 
the intergalactic medium, associated with structure along line 


of sight 12), or the environment of the source CMl- If the 
dispersion is dominated by electrons along the line of sight, 
then the average dispersion-distance relation can be modeled 
1 13l4l5n and the sources appear to be at cosmological dis¬ 
tances of order gigaparsecs. Such cosmological distances per¬ 
mit the study of the clustering of FRB sources in three dimen¬ 
sions, with an expectation that these events trace the large- 
scale structure of the Universe on linear scales. We note that 
while our discussion here is framed in terms of FRBs, other 
bright radio transients that fluctuate on time scales shorter than 
the typical differential lag of ~ Is/100 MHz might also be 
useful as standard pings. 


Neither redshift nor dispersion are perfect proxies for radial 
distance. Redshifts are systematically biased by the peculiar 
velocities of objects relative to the Hubble flow. This leads 
to additional apparent clustering in redshift space which was 
first calculated by Kaiser CS using linear theory. Similarly, 
as we show here, distance estimates from dispersion will be 
biased by inhomogeneities in the electron density Il2l . again 
leading to additional apparent clustering of sources in disper¬ 
sion space. In this Letter we calculate these dispersion-space 
distortions and consider the detectability of the signal by up¬ 
coming surveys. 


Clustering in dispersion space .—The dispersion measure 
of a signal observed in some angular direction h and originat¬ 
ing from comoving radial distance x is 


BM{n,x)=[ dx'a{xfne{nx',x!)- (1) 

Jo 

Here ne{x,x) i^ free electron density as a function 
of location and conformal time. Note that we use x as 
our radial distance and time coordinate, as opposed to red¬ 
shift. We model the cosmological electron density as con¬ 
taining a homogeneous part and perturbations, ne{x,x) = 
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Ueix) [1 + Se{x, x)]- The dispersion measure is thus 
rx 

DM{n,x)= dx'a{x')'^ne{x')[^ + ^e{nx',x')]- (2) 
Jo 

Dispersion space is the three-dimensional coordinates, Xg, 
inferred from the dispersion measure assuming that the elec¬ 
trons are homogeneous. This only affects the radial coordinate 
such that Xg = fiXs, with Xs defined by the equation 


DM(xs) = / dx'aix'frieix')- 


(3) 


Combining the above equations (keeping only terms first 
order in the density perturbations), we find that 


= l + 6e{nx,x), 
dx 


and thus 


(4) 


- X = 


dxSeinx)- 


(5) 


We wish to relate the density of a tracer, /, measured 
in dispersion space to its density in real space. We follow 
the derivation in Kaiser S of the redshift-space distortions. 
Start by noting that the total number of tracers in a volume 
element is the same in both spaces; 


nfg{xg) d^Xg =nf{x)(Px. 


(6) 


We split the density into a homogeneous part plus perturba¬ 
tions, 

n/siXs) [1 + Sfg{xg)] d^Xg = nf{x) [1 + <5/(^)] d^x. (7) 

Averaged over the sky, (xs) = X> ^nd thus the background 
density should be the same in dispersion space as in real space. 


nfsix) = nf{x)- 


Therefore, 


nfsiXs) = nf(x) + {Xs -x) 


dnf 

dx 


= nfix) + 


drif 

dx 


dx'deinx'). 


The Jacobian in spherical coordinates is 




cDxg 


dx X^ 
dXs xl 

l-6e-- [ dx'Seinx')- 
X Jo 


( 8 ) 
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Substituting Egs. fTOl and fT2] into Eq.|7] we obtain 


r r c- , 1 dUf 2 

dfs — df be ( ^ 1 I I I 

^ nf dx xj Jo 


dx'deinx')- (13) 


In the above equation, the —be term is most analogous to the 
Kaiser redshift-space distortions. It is a dilution of tracers in 
dispersion space due to an excess of electrons between the 
tracers. However we note that this term is isotropic in contrast 
to the Kaiser term. This is because any wave vector electron 
perturbation causes dispersion-space distortions, whereas the 
radial velocities that cause redshift-space distortions are only 
sourced by perturbations with radial wave vectors. 

The -J- -JJJ- term arises because the misinterpretation of the 
radial distance causes the observed tracer density to be com¬ 
pared to the wrong background density. The ^ term is caused 
by a misinterpretation of angular distances when the radial 
distance is mismeasured. In both cases, analogous terms are, 
in principle, present in redshift space but are negligible. Be¬ 
cause radial velocities are only sourced by modes with a large 
radial wave vector, there is near perfect cancellation along the 
line of sight, and thus there is very little net error in the radial 
distance. 

Eor brevity in the following sections, we define the coeffi¬ 
cient of the integral term as 


^(x) = 


1 dnf 
Uf dx 


2 

x’ 


(14) 


Large-scale structure is usually studied through its two- 
point statistics, most commonly the power spectrum, P{k). 
Unlike the redshift-space distortions, Eq. [13] does not have a 
simple form in harmonic space. The equation’s third term 
couples harmonic modes, and thus the two-point statistics 
cannot be phrased as a simple power spectrum. We will in¬ 
stead use C'|'*(x,xO’ which is the cross-correlation angular 
power spectrum of the dispersion-space overdensity, on shells 
at X and x': 

Cr(x,x')= / dndn'YUn)YLi.n'){5fs{nx)5fg{n'x:)). 

(15) 


The first two terms in Eq. [13] are stationary. Eor station¬ 
ary, isotropic tracers x and y, we have {Sx{k, x)dy{k' , x)) = 
(27r)^b^(fc — k')Pxy{k, x)- If, for a moment, we ignore struc¬ 
ture evolution, the angular cross-power spectrum of such trac¬ 
ers is 

2 /■“ 

Cr{x,x') = - dkk'^je{kx)je{kx')Pxy{k). (16) 

Jo 

In reality the power spectrum evolves on the order of a Hubble 
time. The angular cross correlations will be very small unless 
X and x' are within a few correlation lengths of one another, 
roughly a hundred megaparsecs. The evolution of the power 
spectrum is negligible over these time differences, which leads 
to a straightforward way to include the evolution; 

2 1 "°° 

cr{x,x') / dkk^ji{kx)je{kx')Pxy{k, (x + x')/2), 
Tl" Jo 

|x-x'l<l/af^- (17) 










The third term in Eq. [T3]is not stationary but is an integral It is straightforward to show that 
over the stationary field Se- Define 5d as 
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^dinx) 



dx'Seinx')- 


(18) 
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0 rx rx' /■°o 

Cfix. x') = - dx” / dx'" / dkejeikx'') 3 dkx’")Pee{k, (x" + x'")/2). 

Jo Jo Jo 

Finally, C|® will contain cross terms between the stationary terms and the integral terms. These will have the form 

o pX 

Ctix, X) = - dx!' / dkk^Ukx')Mkx'')Pe.{k, ix' + x")/2). 

Jo Jo 

Assembling all of these expressions with the proper coefficients, we have 
2 f°° 

Crix, X) = - dkk'^jp{kx)jt{kx')P[ff+ee- 2 ef](k, (x + x')/2) 

Jo 

p rX rX roc 

+ -A{x)A{x') dx" dx'" dkk^Mkx")je{kx"')Pee{k,{x" + x"')m 
^ Jo Jo Jo 

2 rx /*oo 

+ -Aix) dx" dkk'^ji{kx')ji{kx")P[ee-fe]{k,{x'+ x")/‘2) 

Jo Jo 

2 fX 

+ -Aix') dx" dkk‘^ji{kx)j£{kx")P[ee-fe]ik,ix + x")/‘J)- 
^ Jo Jo 


(19) 


( 20 ) 


( 21 ) 


Here, expressions such as P[yj_|_ee- 2 e/] shorthand for 

Pff P Pee — ‘J Pef- 

EquationlSTIcan be simplified substanti ally by adopting the 
small angle and Limber approximations 11 1 8142011 . The small 
angle approximation eliminates the k integral over spherical 
Bessel functions, replacing it with a Fourier transform, and 
is valid for ^ ^ The Limber approximation assumes that 
only modes with a small radial component of their wave vec¬ 
tor contribute to the radial integrals and is valid if the power 
spectra evolve slowly compared to the correlation length Kill 
(which has already been assumed). With these approxima¬ 
tions, we have 


Cfix.x')^ 

^Xmin 1 

+ A{x)A{x') / dx"—Pee{>ylx",x") 

Jo X 

. ^(xmax) r> / / \ 

2 [ee—e/] (^/Xmini Xmin); (22) 

^min 

where = £ + 1/2, Xmin = min(x,x')^ and Xmax = 
max(x,xO ■ We have found that these approximations are 
accurate to within 3% at £ > 10 and use this form for the re¬ 
mainder of the Letter. In the above equation, we dub the three 
terms the “local”, “integral” and “cross” terms, respectively, 
and we will refer to them as such henceforth. 

Modelling and measurement. —The second term in Eq. |22] 


is an integral that mixes many spatial scales: small scales 
from smaller radial distance and larger scales from greater 
distances. The mixing of scales complicates the interpretation 
of measurements since structure formation on scales smaller 
than ~ 10 Mpc/ft. is nonlinear and hard to model. A similar 
issue exists in the field of weak gravitational lensing, where 
the formula for the shear angular power spectrum has a similar 
form. In lensing, a kernel in the line-of-sight integral arising 
from geometric effects suppresses the contribution from small 
distances, partially alleviating the issue. 

However, the lack of a kernel simplifies tomography—the 
use of sources at multiple redshifts to unmix spatial scales. 
Lensing tomography is inexact because the kernel’s shape de¬ 
pends on the source redshift, thus reweighting the line-of-sight 
integral. In contrast, the integrand in the dispersion-space 
integral term is completely independent of the limit of in¬ 
tegration, meaning contributions from different parts of the 
line of sight (and thus different spatial scales) can be sepa¬ 
rated exactly. We leave a study of the potential of dispersion- 
space tomography to future work. To avoid our results be¬ 
ing sensitive to the hard-to-model small scales, we assume 
that the contribution to the integral term from radial distances 
x" < 500 Mpc//i are well measured by the correlations with 
these radial slices. We thus ignore contributions to the inte¬ 
gral term from below this distance and neglect the information 
from these slices in our sensitivity measurements. The small¬ 
est spatial scale that contributes to our plotted angular power 
spectra out to £ = 1000 is thus k = 2 h/Mpc. 
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FIG. 1. (Top) Comoving density of sources in toy model, as de¬ 
scribed in text. Normalization is such that the survey has a total of 
10 000 FRB events. (Bottom) Resulting coefficient A{x) as given in 
Eq.[l4] 


For plotting the dispersion-space power spectra, we use a 
toy model for the electron and FRB clustering that assumes 
both are biased tracers of the dark matter: i.e. Pxy{k) = 
bxbyP{k). We calculate the time dependant dark-matter 
power spectrum using CAMB, using the integrated HALOFIT 
to approximate the nonlinear evolution. We assume that the 
electron bias is he = which should be true on large scales 
due to the approximate conservation of free electrons. The 
bias of FRB sources is unconstrained by data and we choose 
it to be 5/ = 1.3, rou ghly the value for star forming galaxies at 
the relevant redshifts 12211 and providing reasonable contrast to 
the electrons. The coefficient A{x) depends on the density of 
detected FRB sources, nf{x)- We use a simple model where 
all events above a fixed flux are detected, the intrinsic lumi¬ 
nosity function has the form from Schechter with index 
— 1, and events at the cutoff of the luminosity function and at 
a radial distance of x = 2350 Mpc/h are at the flux detec¬ 
tion threshold. This source density is shown in Fig. [T] with the 
resulting dispersion-space power spectrum in Fig. |2l 

The amplitudes of the local and cross terms depend on the 
difference between the electron bias and the highly uncertain 
FRB bias, and the sign of the cross term depends on which of 
these is larger. However, since the integral term dominates for 
most radial distance pairs, our final sensitivity estimates are 
largely unaffected by this uncertainty (as shown below). In 
a high-precision survey, the sensitivity of the cross and local 
terms to the difference in the biases may help break measure¬ 
ment degeneracies, while sign changes of terms as a function 
of radial distance are a signature of bias evolution. 

In a real survey, measurement of the dispersion-space clus¬ 
tering will be complicated by extra contributions to the dis¬ 
persion. First, all signals must pass through the Milky Way’s 
galactic disk, halo, and local environment. However, the re¬ 
sulting dispersion is only a single function of n which pre¬ 
sumably can be well measured. As such, it should be possible 



FIG. 2. Terms in the dispersion-space cross-correlation angular 
power spectrum. All terms are evaluated at x = 2000Mpc//i and 
X — = 10 Mpc//i, except the lowest most local term curve, which 

is for X~ X =50 Mpc//i. While the other two terms are highly 
insensitive to the separation of the radial slices, the local term drops 
rapidly with separation. 


to subtract the local contribution, as was done in Dolag et al. 
i^ . and this is unlikely to be a limiting obstacle. 

Perhaps more concerning is the dispersion from the 
source’s immediate environment as well as the halo in which it 
resides. These electrons will presumably be clustered around 
the source on scales far smaller than the survey resolution and 
are thus not accounted for in the 5^-5 / correlations included 
in the above formalism. This will have a mean contribution, 
which is common to all sources and which may be epoch de¬ 
pendant, that will modify the mean dispersion-distance rela¬ 
tion. This will need to be calibrated in some way or else fit 
with nuisance parameters during parameter estimation. There 
will also be a stochastic piece which varies from source to 
source. This will limit the precision at which Xs can be mea¬ 
sured, thus limiting the resolution of dispersion-space density 
maps in the radial direction. We note that one-point statistics 
from the same data set should provide empirical information 
about these properties independent of their physical origin. 

Detailed treatment of the above source modelling is be¬ 
yond the scope of this Letter. The magnitude of the con¬ 
tribution from the progenitor is, at present, unknown. To 
crudely deal with it here, we use relatively large bin widths 
of Ax = lOOMpc/ft. and exclude the x = xl bin pairs. Be¬ 
cause the local term drops rapidly with radial separations, this 
effectively discards the contribution from the local term which 
is most sensitive to having precise radial distances. 

In first measurements of the angular power spectra, uncer¬ 
tainty will be dominated by the finite number of observed 
events, or shot noise. The noise power spectrum on shells 
at a distance xi with width Ax is 

pN _ 

~ nfixi)x^^x' 


( 23 ) 
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FIG. 3. Sensitivity of a survey with 10 000 dispersion measures 
distributed uniformly over half of the sky. Plotted is the cross¬ 
correlation power spectrum weighted averaged over all pairs of ra¬ 
dial bins. Weights are chosen to maximize signal to noise at f = 100 
for the bf — 1.3 case. 


Here, to indicate that we are now working with quantities 
binned in the radial direction, we have switched to Ce^ij in¬ 
stead of Cg^XjX')- Th® uncertainty on the angular cross¬ 
power spectrum is then 




1 


£(£+l)A£/sky 


, (24) 


where the second term in brackets only contributes for i = j. 

We consider the sensitivity of a survey with a total of 10 000 
dispersion measures from FRB events distributed uniformly 
over half of the sky. To give an idea of the overall sensitivity 
of the survey, we collapse CfX to a single function of i by 
taking a weighted average of all pairs of radial bins over the 
range 500Mpc/h to 3500 Mpc//i. For weights we use the 
signal over noise squared at £ = 100, which maximizes the 
signal to noise at that i and is near optimal at other multipoles. 
To determine if such a survey would be able to distinguish 
between models of FRB clustering we plot the same quantity 
for bf = 0.8. The result is shown in Fig. |3] 

Discussion and conclusions .—We have proposed using 
cosmological FRB sources as standard pings to trace the 3D 
large-scale structure of the Universe using dispersion distance. 
It can be seen in Fig. |2]that the signal is dominated by inho¬ 
mogeneities in the electron density along the different lines of 
sight, inducing misestimates of radial distance, and thus ap¬ 
parent clustering. Because it is the line-of-sight integrated 
structure that correlates, the signal is fairly insensitive to ra¬ 
dial distance and the radial separation of slices, with these de¬ 
pendencies being dominated by the coefficient A[x)- These 
dispersion-space distortions could be used to study the distri¬ 
bution of free electrons in the Universe on large scales. There 
have been other proposals to use FRBs for cosmology 0171124 - 
271] : however these schemes have either not been three di¬ 
mensional or required externally measured redshifts for the 
sources. 

We have shown in Fig.|3]that a survey detecting 10 000 fast 
radio bursts of cosmological origin could detect the cluster¬ 
ing signal. This is in rough approximation to what might be 


achieved in a reasonable amount of time by the Canadian Hy¬ 
drogen Intensity Mapping Experiment (CHIMEfl the Square 
Kilometre ArrajH, or other upcoming wide-held telescopes. 
Such a survey would be primarily sensitive to the integral term 
and thus rather insensitive to the details of ERB clustering 
models. Looking forward, dispersion-space clustering could 
become a precision probe of large-scale structure should ob¬ 
servational factors turn out to be favorable compared to spec¬ 
troscopic surveys. The baryon acoustic oscillation feature can 
be clearly seen in Eig. |2]in the cross term and both local-term 
curves. It may also be possible to extract the baryon acoustic 
oscillation feature from the integral term using tomography. 
This, however, requires that the statistics of the contribution 
to dispersion from the source’s environment and host halo to 
be well understood or small compared to the contribution from 
the line of sight. 

Einally, we note that alone—or in combination with redshift 
data on a subset of the population, independent redshift sur¬ 
veys, or gravitational lensing surveys—dispersion-space sur¬ 
veys could yield very powerful probes of cosmology. Cosmo¬ 
logical dispersion-space data could open a new window into 
the structure of the Universe. 
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